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ABSTRACT 

A  unique  analysis  is  applied  to  the  normal  modes  of  the  U.S.  Navy's  global  operational  numerical  weather 
prediction  model  to  investigate  time-varying  responses  in  the  tropics  and  extratropics  to  tropical  Pacific  sea 
surface  temperature  anomalies.  With  this  new  analysis,  the  modes  are  partitioned  according  to  their  latitudinal 
variances.  This  allows  the  response  energy  to  be  separated  into  tropical  and  extratropical  contributions.  The 
partitioned  responses  are  derived  by  grouping  those  modes  whose  fractional  variance  within  a  present icd  latitudinal 
band  bp  exceeds  some  threshold  value  J.  Since  the  parameters  bp  and  ,i  may  be  chosen  arbitrarily .  this  technique 
gi^atly  increases  the  flexibility  of  the  normal-mode  diagnostic  approach. 

The  partitioned  responses  reveal  distinct  differences  between  the  evolution  and  vertical  scales  of  the  dominant 
modes  in  the  tropics  and  extratropics.  In  the  tropics,  the  structure  is  dominated  by  the  external  mode  and  a 
medium-depth  internal  mode.  The  internal  mode  is  determined  by  the  profile  of  the  large-scale  divergence  and 
subsequent  rotational  wind  ( Walker  circulation )  response  driven  by  enhanced  convection.  In  the  extralropics. 
the  dominant  structure  is  equivalent  barotropie.  The  external  rotational  modes  grow  rapidly  within  the  extralropics 
in  a  manner  that  suggests  that  meridional  propagation  alone  does  not  fully  explain  the  growth  of  the  extratropical 
response. 


1.  Introduction 

This  is  Part  II  of  a  study  investigating  global-scale 
interactions  between  the  tropics  and  extratropics  using 
the  normal  modes  of  a  global  numerical  weather  pre¬ 
diction  (NWP)  model.  The  main  objectives  of  this 
study  are  to  demonstrate  the  rapid  response  of  the 
global  circulation  to  tropical  forcing  anomalies  and  to 
investigate  the  evolution  and  dynamic  structure  of  this 
response  in  a  sophisticated  NWP  model.  The  model 
used  is  part  of  the  Navy  Operational  Global  Atmo¬ 
spheric  Prediction  System  (NOGAPS)  described  by 
Hogan  and  Rosmond  ( 1991 ). 

The  analysis  in  Part  I  (Gelaro  1992,  hereafter  re¬ 
ferred  to  as  GI)  was  conducted  from  a  global  perspec¬ 
tive,  in  which  the  normal  modes,  in  conjunction  with 
standard  difference  field  diagnostics,  were  used  to  an¬ 
alyze  the  model  response  to  sea  surface  temperature 
(SST)  anomalies  in  the  tropical  Pacific.  It  was  shown 
that  the  response  was  dominated  by  external  and  me¬ 
dium-depth  internal  rotational  modes  that  grew  rap¬ 
idly,  obtaining  a  significant  fraction  of  their  steady- 
state  amplitudes  during  the  first  two  weeks  of  the 
simulations.  A  wavenumber  decomposition  of  the 
extratropical  wave  pattern  supported  these  results,  re- 
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vealing  the  rapid  growth  of  the  pattern  between  days 
7  and  10  of  the  simulations.  The  results  demonstrated 
that  the  modes  could  be  used  effectively  to  identify  and 
monitor  the  growth  of  the  dominant  vertical  and  hor¬ 
izontal  structures  of  the  global  response,  and  in  some 
cases,  to  make  conclusions  about  the  physical  processes 
involved.  These  conclusions,  however,  were  limited  by 
the  global  nature  of  the  calculations  presented  in  GI. 

In  this  paper,  the  normal-mode  analysis  is  extended 
to  investigate  separately  the  tropical  and  extratropical 
components  of  the  responses.  This  spatial  segregation 
of  the  response  provides  further  insights  into  the  re¬ 
sponse  mechanisms.  In  sections  2  and  3.  a  new  tech¬ 
nique  is  developed  whereby  the  model  normal  modes 
are  partitioned  according  to  their  latitudinal  variances 
in  order  to  define  the  tropical  and  extratropical  con¬ 
tributions  to  the  response  energy.  The  development  of 
this  technique  follows  directly  from  the  normal-mode 
formulation  of  a  simple  expression  fot  the  total  energy 
in  the  linearized  model,  as  described  in  the  Appendix. 
In  section  4.  this  partitioning  technique  is  used  to  an¬ 
alyze  the  results  of  the  model  experiments  conducted 
in  GI.  The  reader  is  referred  to  sections  2  and  3  of  that 
paper  for  a  description  of  those  experiments  and  an 
overview  of  the  model  response.  In  section  5.  selected 
response  fields  are  projected  onto  the  modes  to  cor¬ 
roborate  the  analysis  and  to  provide  physical  interpre¬ 
tations  of  the  results,  f  inally,  conclusions  and  sugges¬ 
tions  for  future  rcsearcn  ate  discussed  in  section  6. 
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2.  Latitudinal  variance 

As  shown  in  the  Appendix,  the  energy  of  a  single 
normal  mode  J  may  he  written  in  the  form 

=  yajctj  \j,  (I) 

where 


A 


B 


(2) 


and  the  limits  of  integration  are  written  generally  for 
the  purposes  of  the  following  discussion.  As  shown 
schematically  in  ( 2 ),  A,  depends  on  two  factors.  Factor 
A  describes  the  fractional  contributions  of  the  various 
forms  of  energy  to  the  total  energy  of  mode  J.  while 
factor  B  describes  the  latitudinal  variance  of  the  mode. 
As  discussed  in  the  Appendix,  the  structures  of  the 
modes  have  been  normalized  such  that  A  j  =  1  when 
Ej  is  equal  to  the  global  total  energy  of  mode  ./.  In 
that  case.  ( 1 )  reduces  to  Eq.  ( 1 )  in  GI. 

If  we  consider  some  component  of  the  total  energy 
of  mode  J.  E'j  Ej.  then  A  j  can  have  any  value  0 
<  A  j  «  1.  which  may  be  interpreted  as  a  weight  equal 
to  the  ratio  Ej/ E,.  Since  it  is  independent  of  the  data. 
A  j  thus  provides  a  convenient  way  to  examine  the 
energetics  of  the  modes  in  further  detail.  For  example, 
from  factor  A  in  (2)  we  see  that  the  fractional  contri¬ 
bution  of  available  potential  energy  to  the  total  energy 
of  mode  J  is  given  by 

A T)=  I  *JM*.  (3) 

H 


In  the  present  study,  factor  B  in  (2)  is  of  primary  in¬ 
terest  as  it  can  be  evaluated  within  some  latitudinal 
range  bp  to  obtain  the  fraction  of  the  total  variance  of 
mode  ./  within  that  range.  This  information,  in  turn, 
can  be  used  to  define  a  measure  of  the  response  energy- 
in  that  range.  Strictly  speaking,  for  arbitrary  values  of 
bp  there  may  be  nonzero  contributions  from  those 
terms  in  B  for  which  n  E  ft'.  It  can  be  shown,  however, 
that  these  terms  do  not  contribute  significantly  to  the 
variances  and  are  thus  neglected  in  the  present  context. 
A  consequence  of  this  simplification  is  that  the  inte¬ 
grand  in  B  becomes  a  -definite  function  that 

has  equal  value  in  the  same  latitudinal  bands  in  each 
hemisphere.  Therefore,  care  should  be  taken  when  se¬ 
lecting  bp  since;  for  example,  no  distinction  can  be 
made  between  the  variance  of  mode  J  in  the  ranges 
5°-10°N  and  5°-10°S. 

We  seek  away  of  using  B  in  (2)  to  determine  whether 
mode  ./  contributes  to  the  tropical  or  extratropicai  re¬ 
sponse  in  the  experiments  described  in  GI.  For  the  T47 


model  truncation  used  in  this  study,  the  functional 
values  of  E’j  are  defined  on  a  Gaussian  transform  grid 
that  has  72  latitudinal  points.  For  our  purposes,  we 
define  the  tropics  as  extending  eight  Gaussian  latitudes 
from  the  equator  in  both  hemispheres  T  his  corre¬ 
sponds  to  the  latitudinal  belt  between  approximately 
1 9°N  and  I9°S.  The  remaining  area  in  both  hemi¬ 
spheres  is  defined  as  the  extratropics.  Accordingly,  we 
define  for  each  mode  the  discrete  variance  estimates 

44 

Ar*’-  I  I  <4) 

and 

a«m»ip»sv  v  [  /'"■'(  p:)}n.  (5) 

n  i  !./  4* 

for  the  tropics  and  extratropics,  respectively,  where  n 
is  the  weight  given  to  the  functional  value  of  f”,'  at 
Gaussian  latitude  p,  ( Abramouitz  and  Stegun  1965  ). 
with  pt  the  southernmost  latitude  and  p-:  the  north¬ 
ernmost  latitude.  The  summation  over  /  in  (4)  thus 
encompasses  an  equatorially  symmetric  latitudinal 
band  with  poleward  boundaries  at  p2.,  and  pAi  in  the 
Southern  and  Northern  hemispheres,  respectively.  The 
summation  over  i  in  (5)  encompasses  two  separate 
latitudinal  bands  bounded  bv  p,  and  u:>  in  the  Southern 
Hemisphere,  and  by  and  p-,-,  in  the  Northern  Hemi¬ 
sphere.  Note  that,  with  these  definitions,  the  tropical 
and  extratropicai  variance  factors  for  mode  ./  are  re- 
lated  to  one  another  by  A./  =  1  -  A./ 

3.  Partitioning  the  modes 

The  values  of  A and  A;'lropl  can  be  used  effec¬ 
tively  as  a  means  of  classifying  mode  ./  as  tropical  or 
extratropicai.  In  this  way,  an  estimate  of  the  tropical 
energy,  for  example,  can  be  obtained  by  computing 
the  quantity 

/-;"r‘,p’=  V  Ej.  (6) 

.If  T 

where  E,  is  the  total  energy  of  mode  ./  as  given  by 
(A6),  but  the  sum  in  (6)  includes  only  those  modes 
for  which  A.'/"’1’1  exceeds  some  threshold  value  ti  for 
tropical  modes.  This  set  of  tropical  modes  is  denoted 
by  T  and  defined  formally  as 

T  =  {./) a';"’p’  =*  d!.  (7) 

By  defining  the  tropical  energy  as  in  ( 6 ),  we  have  used 
the  values  of  Al.'mp)  as  selection  criteria  for  tropical 
modes  rather  than  as  weighting  functions  a  sum¬ 
mation  over  all  modes  as  might  be  inferred;  for  ex¬ 
ample,  from  ( I  ).  In  this  way.  a  shar|>cr.  but  still  real¬ 
istic,  separation  is  obtained  between  the  characteristics 
of  the  tropical  anu  extratropicai  eoiuubulions  to  the 
global  energy.  Also,  this  approach  is  somewhat  easier 
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FiC.  1.  Values  of  the  variance  factors  A  /  for  the  external  rotational  modes  plotted  as  a 
function  of  their  natural  periods.  The  values  along  the  ordinate  correspond  to  the  fraction  of  the 
modal  variance  between  l^N  and  19°S. 


to  interpret  physically  since  each  term  in  (6)  accounts 
for  the  total  variance  of  mode  J. 

Clearly,  the  key  to  partitioning  the  modes  as  in  (7) 
is  to  define  a  reasonable  threshold  value  &  for 
identifying  tropical  (or  extratropical )  modes.  Unfor¬ 
tunately,  there  is  no  a  priori  way  of  knowing  what 
threshold  values  of  AyIrop>  and  A  yx  rop)  adequately  sep¬ 
arate  the  two  groups  of  modes.  The  most  straightfor¬ 


ward  approach  is  simply  to  plot  the  variance  factors  as 
a  function  of  the  natural  frequencies  wy  of  the  modes 
in  order  to  see  where  a  natural  separation  might  exist. 

Figures  1  and  2  show  the  values  of  Aylropl  corre¬ 
sponding  to  the  rotational  and  eastward  gravitational 
modes,  respectively,  for  the  external  (/  =  1 )  vertical 
mode.  Again,  these  values  are  computed  for  bn  ex¬ 
tending  from  approximately  19°N  to  19°S,  as  given 
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by  (4).  and  describe  the  fraction  of  the  total  variance 
of  mode  J  that  occurs  within  this  latitudinal  band.  Ihc 
frequencies  are  ordered  from  smallest  to  largest  ac¬ 
cording  to  their  absolute  values.  This  ordering  is  from 
largest  to  smallest  spatial  scales  for  the  grav  itational 
modes,  and  from  smallest  to  largest  spatial  scales  for 
the  rotational  modes.  As  in  Gl.  we  consider  all  modes 
except  the  zonally  symmetric  ones  corresponding  to  m 
-  0.  In  both  figures,  there  appears  to  be  a  lower  cutoff 
near  A'/>P  -  0---  indicating  that  very  few  modes  have 
less  than  20r<  of  their  variance  in  the  tropics.  In 
other  words,  there  are  no  purely  extratropical  modes 
[i.e..  values  at,  or  near.  A“ri,p)  -  0]  according  to  this 
definition.  In  contrast,  note  that  in  both  figures  there 
are  modes  for  which  Aylrop)  is  very  close  to  one,  indi¬ 
cating  that  nearly  all  of  their  variance  occurs  within 
the  tropics.  The  value  of  the  lower  cutoff  depends,  in 
part,  on  the  choice  of  ?>n  and  increases  or  decreases, 
accordingly.  In  both  figures,  the  values  of  Ay rop)  appear 
to  be  organized  into  discrete  bands  that  arc  downward 
from  left  to  right  for  the  rotational  modes  in  Fig.  1, 
and  downward  from  right  to  left  for  the  eastward  grav¬ 
itational  modes  in  Fig.  2.  Each  band  represents  a  family 
of  modes  having  the  same  meridional  scale  denoted 
by  the  index  j,  while  each  mode  in  that  band,  or  j 
family,  corresponds  to  a  different  zonal  wavenumber 
m.  There  are  47  such  families  (excluding  the  zonally 
symmetric  modes)  corresponding  to  the  T47  model 
truncation,  in  which  the  number  of  modes  varies  from 
47  for  7=1,  down  to  l  for  j  =  47. 

For  the  rotational  modes  in  Fig.  1.  the  7  =  1  family 
is  the  upper-rightmost  band  and  represents  the  mixed 
Rossby-gravity  modes.  Many  of  these  modes  are  trop¬ 
ically  trapped  as  indicated  by  their  relatively  large 


values  of  Ay""1' ”.  Ihc  /  1  family  for  the  eastward 

gravitational  modes  is  the  upper-leftmost  band  in  Fig. 
2  and  corresponds  to  low-frequency  Kelvin  modes  As 
expected,  these  modes  also  tend  to  be  trapped  and  have 
values  close  to  A;"'1"  I.  Note  that  in  both  figures 
those  modes  with  large  meridional  structure  (small 
values  of/)  have  more  of  their  variance  in  the  t.opics 
than  those  with  smaller  meridional  structure  ( large 
values  of  /).  Clearly,  there  is  much  more  information 
that  can  be  gleaned  from  plots  of  A ,  for  various  ranges 
of  bp.  For  example,  it  is  interesting  to  note  that  the 
rotational  modes  in  Fig.  !  exhibit  very  little  tropical 
variance  for  wy  <  2.3  X  10  "  s  1  ( periods  greater  than 
31  days).  At  this  time,  it  is  unclear  whether  there  is 
some  theoretical  basis  for  this  sharp  cutoff. 

Plots  such  as  those  in  Figs.  1  and  2  for  the  remaining 
vertical  modes  reveal  similar  patterns  for  their  corre¬ 
sponding  frequency  ranges.  As  an  example,  we  show 
in  Fig.  3  values  of  Ay""1”  corresponding  to  the  eastward 
gravitational  modes  for  /  =  5.  For  all  vertical  modes, 
it  was  generally  found  that  the  modes  demonstrated  a 
marked  increase  in  tropical  variance  for  /  •?  6.  Roughly 
speaking,  most  of  the  modes  in  these  families  have 
values  of  >  0.4.  Based  on  these  figures,  we 

adopted  the  cutoff  value  A yUl,p;  5  0.4  for  "tropical” 
modes,  and  Ay"01”  ?  0.3  for  “extratropical"  modes. 
Those  modes  for  which  0.3  <  A <  0.4  are  difficult 
to  classify  as  belonging  in  one  group  or  the  other,  and 
are  thus  excluded  from  the  analysis.  Note  that  it  is  only 
necessary  to  consider  Ay"°pl  in  order  to  identify  both 
the  tropical  and  extratropical  modes  since  Ay'"’p'  and 
Ay"rop>  sum  to  one.  as  discussed  earlier:  that  is.  Ay""1” 
<  0.3  is  synonymous  with  Aj'mpl  >  0.7. 

It  is  important  to  note  that,  for  the  threshold  values 
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Fig  .  3.  As  in  Fig.  I  except  for  the  1-5  eastward  gravitational  modes. 
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Table  I.  The  numbers  of  rotational  (RT).  eastward  gravitational 
(EG),  and  westward  gravitational  |WG)  modes  classified  as  "tropical" 
or  "extratropical"  for  each  vertical  mode  /  based  on  the  threshold 
values  A'}”*'  >  0.4  and  Aj',r'pi  >  0.7.  respective!). 
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Tropical 

Extratropical 

RT 

EG 

WG 

RT 

EG 

WG 
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748 
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245 
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612 

14 

254 

248 

253 

630 

643 
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16 
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263 

630 
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598 

17 

25? 

256 

265 

632 

634 

596 

IS 

250 

255 

271 

63.3 

634 

59.3 

calculations  shown  in  Cil.  we  plot  for  each  vertical 
inode  a  normalized  response  energy  that,  tor  the  trop¬ 
ical  energy  in  Fig.  4a.  takes  the  foim 


where  the  summation  includes  only  tropical  mixles  that 
have  equivalent  depth  /(/(denoted  by  the  set  T,)  and. 
for  notational  convenience,  the  superscript  b  used  in 
Gl  to  denote  a  response  held  is  implied.  I  or  the  extra- 
tiopieal  energy,  we  plot  tor  each  vertical  mode  the 
quantity  w  hich  we  define  by  analogy  to  ( 8  I  by 

summing  over  only  the  set  of  extratropical  modes  Jl 
and  normalizing  by  the  extratropical  energy  r  . 

The  proliles  in  both  regions  are  qualitativ  ely  similar 
during  days  1-5.  w  ith  relative  maximums  at  /  4  and 

/  =  16.  For  the  extratropical  modes  in  Fig.  4b.  the  ex¬ 
terna)  modes  are  also  relatively  energetic  during  this 
time.  By  days  21-25.  the  external  modes  dominate  the 
profiles  in  both  regions,  but  the  relative  contributions 


selected  above,  the  number  of  horizontal  modes  of  each 
type  is  fairly  evenly  distributed  within  each  group.  Ta¬ 
ble  1  shows  the  numbers  of  rotational  and  eastward 
and  westward  gravitational  modes  in  the  tropical  and 
extratropical  gtoups  for  each  vertical  mode.  Roughly 
speaking,  we  see  that  for  /  <  8  there  are  between  200 
and  220  tropical  modes  of  each  type,  and  between  710 
and  740  extratropical  modes  of  each  type  for  each  ver¬ 
tical  mode.  This  ensures  that  a  fair  comparison  can  be 
made  between  the  relative  contributions  of  each  vertical 
mode  to  the  tropical  or  extratropical  response  energy. 
There  are.  of  course,  more  extratropical  modes  than 
tropical  modes  because  of  the  difference  in  the  latitu¬ 
dinal  ranges  bn  used  to  define  each  region.  For  each 
vertical  mode,  the  two  groups  account  for  roughly  85G 
of  the  total  number  of  modes  in  the  complete  ( global ) 
set.  That  is.  15r4  of  the  complete  set  have  values  0.3 
<  Ayr"pl  <  0.4.  and  are  not  included  in  either  group. 
Further  justification  for  the  threshold  values  of 
A,/lmpl  and  A 'j'r"p)  selected  above  are  given  in  the  fol¬ 
lowing  sections. 

4.  Tropical  and  extratropical  responses 

In  this  section,  the  normal-mode  partitioning  tech¬ 
nique  developed  in  sections  2  and  3  is  used  to  analyze 
the  tropical  and  extratropical  contributions  to  the 
global  response  energy  in  the  western  Pacific  SSI 
anomaly  experiment  described  in  GI.  Figure  4  shows 
the  response  energy  in  each  region  for  all  mode  types 
as  a  function  of  vertical  mode  /  tor  the  two  5-da\ 
mean  periods  corresponding  to  days  I  -5  and  days  2  I  - 
25  of  the  simulations.  Analogous  to  the  global  energy 
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from  the  medium-depth  internal  modes  now  differ 
markedly.  In  the  tropics,  /  =  4  remains  the  dominant 
internal  mode,  while  /  =  5  now  dominates  in  the  ex¬ 
tratropics.  In  both  regions,  only  modes  corresponding 
to  /  <  6  contribute  significantly  by  this  time. 

The  evolution  of  the  rotational-mode  energy  in  the 
tropics  and  extratropics  for  vertical  modes  1-5  is  shown 
in  Fig.  5.  Again,  the  external  modes  dominate  in  both 
regions,  but  the  importance  of  /  •=  4  in  the  tropics  and 
/  =  5  in  the  extratropics  is  clearly  evident.  Interestingly, 
the  rotational  modes  appear  to  grow  more  rapidly  in 
the  extratropics.  This  result,  combined  with  the  results 
in  Fig.  1 2  in  GI,  w  hich  showed  that  certain  modes  grow 
at  preferred  locations  in  the  midlatitudes,  suggests  that 
meridionally  propagating  tropical  wave  energy  does  not 
completely  explain  the  extratropical  response.  Also, 
note  the  relatively  smooth  evolution  of  the  dominant 
extratropical  modes  compared  with  those  in  the  tropics. 
Figure  6  shows  that  the  gravitational  modes  play  a  more 
important  role  in  the  overall  energy  balance  in  the 
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Fig.  6.  As  in  Fig.  5  except  for  the  gravitational  modes.  Note  the 
difFerent  ordinate  scale  and  smaller  energies  than  in  Fig  5. 


tropics  but,  in  the  long  term,  arc  still  roughly  an  order 
of  magnitude  smaller  than  the  rotational  modes. 

Since,  strictly  speaking,  the  results  in  Figs.  4-6  reveal 
only  information  about  the  structure  of  the  responses, 
an  interesting  question  is  w  hethcr  the  predominant  re¬ 
sponse  modes  in  the  tropics  and  extratropics  arc  ac¬ 
tually  of  different  dynamical  origins.  The  predomi¬ 
nance  of  the  external  modes  in  Fig.  5b  is  consistent 
with  the  equivalent  barotropic  structure  of  the  extra¬ 
tropical  response  (cf.  Fig.  8  in  GI ).  In  a  strongly  sheared 
environment,  this  response  has  its  maximum  ampli¬ 
tude  near  tropopause  level  (e  g..  Held  1983).  As  dis¬ 
cussed  in  GI,  however,  and  shown  here  in  Fig.  7.  the 
linearization  used  in  the  present  study  yields  an  external 
mode  that  is  nearly  barotropic  throughout  the  tropo¬ 
sphere  and  thus  is  inadequate  for  describing  this  struc¬ 
ture  completely.  Thus,  the  simultaneous  evolution  of 
/  =  5  in  the  extratropics,  which  has  its  maximum  am¬ 
plitude  at  approximately  250  mb  ( Fig.  7  >.  undoubtedly 
is  also  a  manifestation  of  the  growth  of  the  equivalent 
barotropic  response.  This  conclusion  is  confirmed  in 
section  6. 
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Fig.  7.  Vertical  modes  I,  4,  and  5  of  the  NOGAPS  model,  based  on 
the  parameter  values  shown  in  Table  I  of  GI  (signs  are  arbitrary). 

In  contrast  with  the  extratropics,  there  is  no  obvious 
barotropic  phenomenon  in  the  tropics  that  explains 
the  strong  external-mode  response  in  Fig.  5a.  A  growth 
mechanism  for  these  modes,  however,  has  been  dem¬ 
onstrated  by  Lim  and  Chang  ( 1986)  (hereafter  referred 
to  as  LC),  whereby  vertical  shear  acts  to  produce  an 
external  response  from  purely  internal  (e.g.,  convec¬ 
tive)  forcing.  LC  argue  that  internal  forcing  initially 
produces  an  internal-mode  response  with  upper-level 
divergence  and  lower-level  convergence  that,  by 
stretching  the  lower  atmospheric  column  and  com¬ 
pressing  the  upper  one,  leads  to  an  internal-mode  ro¬ 
tational  response  through  conservation  of  potential 
vorticity.  Over  a  period  of  several  days,  vertical  wind 
shear  acts  to  couple  the  internal  and  external  rotational 
modes,  resulting  in  a  transfer  of  energy  to  the  latter. 
With  this  in  mind,  it  is  interesting  to  note  that  the  / 
=  4  rotational  response  in  Fig.  5a  is  preceded  by  a 
strong  1  =  4  gravitational  (i.e.,  divergent)  response,  as 
shown  in  Fig.  6a.  The  projection  of  the  gravitational 
response  onto  /  —  4  is  explained  by  Fig.  8a,  which  shows 
a  vertical  profile  of  the  30-dav  mean  divergence  re¬ 
sponse  averaged  over  the  region  of  anomalous  SST  in 
the  GI  experiments.  The  structure  of  the  profile,  char¬ 
acterized  by  a  sharp  peak  with  maximum  divergence 
just  below  100  mb,  and  a  deep,  but  weaker,  layer  of 
convergence  below  200  mb  that  extends  to  the  surface, 
is  very  similar  to  that  of  vertical  mode  4  (see  Fig.  7; 
signs  are  arbitrary).  As  suggested  by  LC,  a  rotational 
response  develops  with  similar  vertical  structure,  as  in¬ 


dicated  by  the  growth  oft  --  4  in  Fig.  5a  and  supported 
by  Fig.  8b,  showing  the  vertical  profile  of  the  30-das 
mean  rotational  zonal-wind  response  averaged  between 
1 1  °N-1 1  °S  and  90°-I70°W.  The  profile  in  Fig.  8b 
reveals  the  enhanced  Walker  circulation  in  the  central 
and  eastern  Pacific  with  strong  westerlies  at  upper  levels 
and  weaker  easterlies  below.  A  response  with  similar 
vertical  structure,  but  opposite  sign,  occurs  to  the  west 
of  the  SST  anomaly  in  the  Indian  Ocean  ( not  shown ). 

Finally,  note  that  in  contrast  with  the  /  =  5  response 
in  Fig.  5b  the  /  =  4  response  in  Fig.  5a  grows  slowly 
and  then  fluctuates  significantly  about  some  quasi¬ 
steady  state.  The  time  scale  of  these  fluctuations  is 
comparable  with  that  observed  in  time  series  of  the 
model’s  convective  rainfall  (Hogan  1991.  personal 
communication).  The  external  modes  in  Fig.  5a  also 
grow  slowly  at  first,  but  then  grow  much  more  rapidly 
after  the  second  week  of  the  simulation  when  the 
growth  at  /  =  4  has  leveled  off.  Despite  the  delayed 
rapid  growth  of  the  external  modes,  it  is  puzzling  that 
there  is  no  observable  lag  in  their  initial  growth  com¬ 
pared  with  1=4.  This  detail  may  be  lost  either  as  a 
result  of  the  partitioning  procedure  [specifically,  the 
choices  for  the  latitudinal  ranges  used  to  define  (4) 
and  ( 5 ),  or  the  threshold  value  0  in  ( 7 )]  or  the  temporal 
resolution  of  these  results.  Another  possibility  is  that 
the  theory  proposed  by  LC  is  somewhat  oversimplified, 
having  been  derived  from  a  much  simpler  model. 
Nonetheless,  the  growth  rate  and  magnitude  of  the  ex¬ 
ternal  response  agree  well  with  those  obtained  by  LC 
and  the  contrasting  character  of  the  growth  curves  in 
Figs.  5a  and  5b  seems  to  indicate  that  the  predominant 
responses  in  each  region  are  of  different  dynamical  or¬ 
igins. 

5.  Vertical-mode  projections 

As  indicated  above,  the  results  in  Figs.  4-6  depend 
on  the  variance  factors  (4)  and  (5),  and  in  particular, 
on  the  threshold  value  0  in  (7)  used  to  define  the  subsets 
of  tropical  and  extratropical  modes.  In  this  section,  the 
upper-tropospheric  streamfunction  response  is  pro¬ 
jected  onto  the  model  normal  modes.  These  projections 
are  done  to  check  the  physical  plausibility  of  the  par¬ 
titioned  responses  discussed  in  section  5.  The  projec¬ 
tions  shown  here  are  totally  independent  of  the  parti¬ 
tioning  technique  used  to  obtain  those  results.  The 
contribution  to  the  response  from  each  vertical  mode 
is  computed  as  the  difference  between  the  total  response 
field  and  a  response  field  for  which  the  normal-mode 
amplitudes  corresponding  to  that  vertical  mode  have 
been  set  to  zero;  that  is,  a  type  of  linear  initialization 
is  performed.  It  was  decided  to  use  this  difference 
method  rather  than  compute  the  desired  projections 
directly  because  of  the  potential  difficulty  in  recon¬ 
structing  the  fields  from  only  a  single  vertical  mode. 
Projections  of  the  streamfunction  response  at  1 50  and 
300  mb  are  presented  since  they  are  close  to  the  levels 
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Fig.  8  Area-averaged  30-day  mean  response  profiles  of  (a)  divergence  in  the  region  of  .ncreased  SST  and  (b) 
rotational  zonal-wind  speed  over  the  central  and  eastern  tropical  Pacific  in  the  G1  experiment  The  profile  in  (b)  may 
be  thought  of  as  representing  the  Pacific  Walker  circulation.  Note  the  similarity  of  both  profiles  to  the  structure  of 
vertical  mode  4  shown  in  Fig.  7. 


of  maximum  amplitude  for  vertical  modes  4  and  5, 
respectively  (see  Fig.  7). 

Figure  9  shows  the  30-day  mean  stu  amfunction  re¬ 
sponse  at  300  mb.  The  total  response  is  shown  in  Fig. 
9a,  while  the  projections  onto  vertical  modes  1, 4,  and 
5  are  shown  in  Figs.  9b,  9c,  and  9d,  respectively.  As 
expected,  the  total  response  in  Fig.  9a  shows  a  strong 
extratropical  stationary  wave  pattern  in  both  hemi¬ 
spheres  and  some  indication  of  anomalous  equatorial 
westerlies  in  the  central  Pacific.  The  extratropical  wave 
pattern,  particularly  in  the  Northern  Hemisphere,  pro¬ 
jects  strongly  onto  vertical  mode  1  (Fig.  9b)  as  well  as 
vertical  mode  5  (Fig.  9d),  in  agreement  with  the  results 
in  Fig.  5b  for  the  extratropical  response  energy.  Fur¬ 
thermore,  the  results  in  Fig.  9  confirm  that  the  growth 
of  /  =  1  and  /  =  5  in  Fig.  5b  is,  in  fact,  a  manifestation 
of  the  equivalent  barotropic  structure  of  this  response. 
For  example,  note  that  most  of  the  variance  in  the 
circulation  anomalies  over  North  America  in  Fig.  9a 
is  accounted  for  by  the  projections  in  Figs.  9b  and  9d. 
In  the  tropics,  there  is  a  weak,  but  clearly  identifiable, 
projection  of  the  equatorial  Pacific  westerlies  in  Fig. 
9b,  while  the  projection  in  Tig.  9d  has  very  little  am¬ 
plitude  anywhere  in  the  tropics.  This  is  consistent  with 
the  relatively  weak  contribution  from  /  =  5  in  Fig.  5a. 
As  expected,  the  projection  onto  /  =  4  in  Fig.  9c  has 


considerably  less  amplitude  anywhere  at  300  mb  in 
accordance  with  the  vertical  structure  of  this  mode  as 
shown  in  Fig.  7. 

Figure  10  is  the  same  as  Fig.  9,  except  for  a  30-day 
mean  streamfunction  response  at  1 50  mb.  Again,  the 
total  response  in  Fig.  10a  shows  a  strong  stationary 
wave  pattern  in  both  hemispheres,  but  also  has  con¬ 
siderably  more  amplitude  in  the  tropics  than  the 
streamfunction  response  at  300  mb.  This  is  consistent 
with  the  strong  response  at  /  =  4  (which  has  its  max¬ 
imum  amplitude  near  this  level)  in  Fig.  5a  and  with 
the  fact  that  the  tropical  isponse  has  its  maximum 
amplitude  near  tropopause  level.  In  Fig.  10a,  anoma¬ 
lous  westerlies  to  the  east  of  the  forcing  region  are 
clearly  defined  by  the  strong  cross-equatona!  st.eam- 
function  gradient  that  extends  across  most  of  the  Pa¬ 
cific,  while  there  are  strong  anticyclones  to  the  north¬ 
west  and  southwest  of  the  forcing.  The  latter  features 
are  well-known  responses  to  an  cquatorially  symmetric 
heat  source  (e.g.,  Blackmom  et  al.  1983:  Palmer  and 
Mansfield  1986:  Sardeshmukh  and  Hoskins  1988).  The 
extratropical  wave  pattern  projects  strongly  onto  ver¬ 
tical  mode  1  in  Fig.  10b,  while  the  tropical  response 
features  have  only  weak  projection  onto  this  mode. 
Note  the  similarity  between  the  character  and  ampli- 
‘ude  of  the  projections  in  Figs.  9b  and  1 0b,  reflecting 
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Fig.  9.  1  he  til- day  mean  300- mh  stream  function  response  in  the  ( II  experiment  lor  i  a )  the  total  rc'pnnsv  and  lor  the  project1  >i.s  of  the 
response  onto  (b)  vertical  mode  I .  (el  vertical  mode  4,  and  Id)  vertical  mode  5  I  he  contour  interval  is  2  5  ■  III'  m  s  with  negative 
values  denoted  by  dashed  lines. 


the  nearly  barotropic  structure  of  the  external  modes.  Again,  these  results  clearly  support  the  results  in  Fig. 

In  contrast.  Fig.  10c  shows  that  the  tropical  response  5a  showing  the  strong  contribution  from  /  4  to  the 

projects  strongly  onto  vertical  mode  4,  while  the  mid-  tropical  response  energy.  Conversely,  the  relative  lack 

die-  and  high-latitude  response  has  much  less  amplitude  of  amplitude  in  tne  extratropics  in  Fig.  10c  is  in  agree- 

for  this  vertical  mode.  Note,  for  example,  the  near  ah-  ment  v>  ith  the  relatively  weak  contribution  from  /  4 

sence  of  the  North  American  response  in  Fig.  10c.  to  the  cxtratropical  response  energy  in  Fig.  5h.  As  ex- 
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peeled,  the  projection  onto  vertical  mode  5  in  Fig.  I<kl 
is  negligible  at  1  ?()  mb.  In  summary,  the  results  in  1  igs. 
9  and  10  support  the  partitioned  responses  in  Fig.  5. 

6.  Summary  and  conclusions 

( ilobal-scale  interactions  between  'he  tropics  ami 
extratropics  acre  examined  in  terms  •>!  the  norma! 


modes  of  a  sophisticated  NWP  model  \  new  technique 
was  developed  whereby  the  modes  were  partitioned 
lecordmg  to  their  latitudinal  variances  in  order  to  de- 
line  the  tropical  and  exlratropica!  components  of  the 
model  total  energy.  I  he  development  of  this  technique 
Followed  From  a  straightforward  derivation  oi  the  model 
total  energy  in  normal-mode  Form,  in  which  the  lati¬ 
tudinal  structures  of  the  modes  appeared  .e  explicit 


Fig.  10.  As  in  Fig.  9  except  at  150  mb. 


terms  in  an  integral  equation.  The  partitioned  responses 
were  derived  from  those  modes  whose  fractional  vari¬ 
ance  within  either  region  exceeded  a  predefined 
threshold  value.  We  applied  this  partitioning  technique 
to  analyze  the  results  of  a  series  of  experiments  con¬ 
ducted  by  Gelaro  ( 1992),  who  used  the  U  S.  Navy’s 
operational  global  forecast  model  to  demonstrate  the 
rapid  atmospheric  response  to  SST  anomalies  in  the 


tropical  Pacific.  The  analysis  revealed  the  different 
structural  characteristics  of  the  tropical  and  extratrop- 
ical  responses  obtained  in  those  expen mcnts.  and  pro¬ 
vided  insight  into  the  mechanisms  that  govern  their 
evolution.  Furthermore,  it  showed  that  the  normal 
modes  are  a  powerful  and  flexible  tool  for  diagnosing 
the  behavior  of  complicated  models. 

Both  external  and  internal  rotational  modes  were 


Fig.  10.  (Continued) 


shown  to  contribute  significantly  to  the  tropical  and 
extratropical  responses.  Of  particular  interest  was  that 
the  dominant  internal-mode  responses  projected  onto 
different  vertical  modes  in  each  region,  with  vertical 
mode  4  dominating  in  the  tropics  and  vertical  mode 
5  dominating  in  the  extratropics.  This  distinction,  in 
addition  to  the  different  temporal  behavior  of  the  trop¬ 
ical  and  extratropical  modes,  was  indicative  of  the  dif¬ 


ferent  dynamical  origins  of  the  dominant  modes  in 
each  region. 

In  the  extratropics,  the  internal  and  external  modes 
evolved  simultaneously  throughout  the  simulations, 
with  the  latter  clearly  dominating  the  response.  Ver¬ 
tical-mode  projections  of  the  upper-tropospheric 
streamfunction  response  confirmed  that  both  modes 
were  required  to  resolve  the  strongly  equivalent  baro- 
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tropic  structure  of  the  extratropical  response,  which 
obtained  its  maximum  amplitude  near  tropopause 
level.  As  in  the  global  analysis  presented  in  Ciclaro 
(  1992  ),  the  response  evolved  rapidly.  with  much  of  its 
long-term  character  well  developed  alter  the  first  week. 
The  projection  of  this  response  onto  vertical  mode  5 
(which  has  its  maximum  amplitude  near  25(1  mb.  hut 
changes  sign  in  the  middle  troposphere )  does  not  itself 
appear  to  be  of  any  dynamical  significance.  Rather, 
this  seems  to  be  an  indication  that  the  no-motion  basic 
state  used  to  derive  the  modes  is  a  poor  approximation 
of  the  model's  very  energetic  and  strongly  sheared 
midlatitude  basic  state.  This  in  no  way  invalidates  the 
applicability  of  normal  modes  deriv  ed  from  simple  ba¬ 
sic  states,  as  demonstrated  by  their  successful  use  in 
nonlinear  normal-mode  initialization  schemes  (c.g.. 
Andersen  1977;  Hogan  el  al.  1991 ).  as  well  as  in  model 
diagnostic  studies  ( e.g..  Errieo  et  al.  1 9KS  ).  It  does  im¬ 
ply.  however,  that  care  must  be  taken  when  interpreting 
the  physical  mechanisms  that  produce  a  particular 
modal  response. 

In  contrast,  the  simple  basic  state  used  in  these  cal¬ 
culations  is  a  considerably  better  ( though  not  optimal  1 
approximation  of  the  model's  tropical  basic  state  Ac¬ 
cordingly.  the  predominance  of  vertical  mode  4  in  the 
tropics  was  shown  to  be  dy  namically  significant,  being 
entirely  consistent  with  the  structure  of  the  large-scale 
divergence  and  zonal-w  ind  fields  there.  The  strong  re¬ 
sponses  at  this  vertical  scale  in  both  the  gravitational 
and  rotational  modes  appeared  to  agree  w  ith  the  results 
of  l.im  and  Chang  (  1986 ).  who  showed  that  convective 
heating  can  produce  a  hierarchy  of  internal-mode  re¬ 
sponses  culminating  in  a  transfer  of  energy  to  the  ex¬ 
ternal  rotational  modes.  However,  it  was  not  possible 
to  verify  this  completely  from  the  results  presented  here. 
Nonetheless,  the  relatively  slower  and  more  variable 
evolution  of  the  rotational  response  in  the  tropics  was 
in  clear  contrast  with  the  smooth  and  surprisingly  rapid 
evolution  of  the  extratropical  response. 

It  is  important  to  bear  in  mind  that  the  results  of 
this  analysis  (as  well  as  the  success  of  the  partitioning 
method  itself)  depend  on  the  latitudinal  ranges  nt u  used 
to  define  the  variance  factors  A',""1’  and  A1,'"’1’1,  and 
on  the  threshold  values  ti  used  to  define  the  tropical 
and  extratropical  modes.  Although  the  values  chosen 
for  these  parameters  yielded  physically  reasonable  re¬ 
sults  (as  verified,  for  example,  by  the  strcamfunction 
projections  in  section  6).  additional  work  is  required 
to  investigate  the  robustness  of  these  results  and  their 
sensitivity  to  the  aforementioned  parameters.  In  par¬ 
ticular.  by  choosing  geographically  separated  latitudinal 
ranges  to  define  the  tropics  and  extratropics  ( rather 
than  having  a  common  latitudinal  boundary  at  19°  in 
each  hemisphere ),  or  by  using  more  stringent  threshold 
values  for  grouping  the  modes,  wc  might  minimize  the 
risk  of  "smearing”  the  responses  in  each  region.  Other 
improvements,  including  additional  simulations  and 


greater  sampling  resolution  than  the  5-day  means  pre¬ 
sented  here,  might  also  resolve  some  of  these  issues 
Despite  these  limitations,  the  results  clearly  demon¬ 
strate  the  usefulness  of  the  normal  modes  as  diagnostic 
tools  that,  because  of  the  inclusion  of  normal-mode 
initialization  schemes  in  most  global  forecast  models, 
can  be  readily  developed. 
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Normal-Mode  Energetics  in  NOCi APS 

In  a  continuous  linear  system,  the  total  kinetic  plus 
available  potential  energy  /.  per  unit  mass  can  be  writ¬ 
ten 


/  l/fjJ \(ln,  (  A  I  ) 

where  is  the  strcamfunction.  X  is  the  velocity  poten¬ 
tial.  <I>  is  the  geopotential,  and  <  is  a  differential  operator 
of  the  vertical  coordinate  n  ( Errieo  1987  ).  The  variables 
A  and  g  are  the  longitude  and  sine  of  the  latitude,  re¬ 
spectively.  and  V'  is  the  horizontal  l.aplaeian  operator 
in  spherical  coordinates. 

For  our  purposes,  we  seek  the  normal-mode  trans¬ 
form  of(  A  I  ).  Because  the  NOCiAPS  model  is  spectral 
in  both  horizontal  dimensions  and  uses  vorticity  C and 
divergence  I)  as  the  dependent  variables  for  the  velocity 
field  (  Hogan  and  Rosmond  1 99 1  ) .  an  appropriate  form 
of  ( A I  )  for  this  study  is 


i - />:»v< Cf:;'*  •  />: .ik;*  )]/*;:' r:  <.\2> 

in  which  C" .  /•>’!  .  and  <I\|'  are  spherical  harmonic  coef¬ 
ficients;  V";  is  an  associated  Legendre  polynomial  of 
degree  m  and  order  n:  r„  :  a'  /  nt  n  t  1  )  derives  from 
the  spectral  form  of  the  l.aplaeian  operator,  where  a 
is  the  radius  of  the  earth;  and  an  asterisk  denotes 
a  complex  conjugate.  The  factors  h  V 2 If  and  < 
V2 52  V  are  dimensional  constants  in  which  52  is  the 
angular  velocity  of  the  earth.  In  obtaining  (  A2)  we 
have  utilized  the  orthogonal  relationship  between 
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Fourier  functions  with  unlike  indices  and  the  condition 
that;  for  example,  =  f"*  for  real  f . 

To  obtain  the  normal-mode  representation  of(  A2), 
we  expand  the  dependent  variables  in  the  form 

( H7,  <fC)T  =  Z  Z  ajSiH'z,,  (  A3 ) 

/  j 


tributions  from  only  those  terms  for  which  n  =  m'  and 
(A5)  reduces  to  the  simple  form 

Ej  =  7 itja* .  ( A6 ) 
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